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We show that the universal tail in the momentum distribution of dilute Fermi gases implies 
that the spectral function A(k,aj) must have weight below the chemical potential for large momen- 
tum k ^ kp, with observable consequences in RF spectroscopy experiments. We find that this 
incoherent spectral weight is centered about uj ~ — e(k) in a range of energies of order vpk. This 
"bending back" in the dispersion, while natural for superfluids, is quite surprising for normal gases. 
This universal structure is present in the hard-sphere gas as well as the Fermi liquid ground state 
of the highly imbalanced, attractive gas near unitarity. We argue that, even in the BCS superfluid, 
this bending back at large k is dominated by interaction effects which do not reflect the pairing gap. 
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The spectral function A(k, oj) = — IniG'(k, w -I- i0+)/7r 
of the single-particle Green's function G is of fundamen- 
tal interest in many-body physics [H, HI • In addition to 
information about the spectrum and dynamics of single- 
particle excitations, it is also directly related to ther- 
modynamic functions of a many-particle system. Very 
recently there has been experimental progress in measur- 
ing (the occupied part of) A(k, uj) in strongly interacting 
Fermi gases 0, 0] , using a momentum- resolved version [5[ 
of radio frequency (RF) spectroscopy These mea- 

surements Q of A(k, w) for ultracold atomic gases are the 
analog of angle-resolved photoemission, which has given 
deep new insights into novel materials. 

In this paper we uncover remarkable universal large- 
k structure of A{k,uj) for dilute gases with observable 
effects in RF experiments. Our investigation was mo- 
tivated by the elucidation of the universal ultraviolet 
structure of equal-time correlations by Tan [13] . One 
of his central results is the universal k ^ kp behavior 
of the momentum distribution n^Ck) ~ C/k^, where 
C is the "contact" H, Using the T = sum 

rule (ia;A(k, w) = n(k), this necessarily implies that 
A(k, uj) has weight below the chemical potential {lu < 0) 
for k ^ kp- This is "incoherent" spectral weight, not 
associated with the coherent Landau quasiparticle. 

We ask the question: Where is this incoherent spec- 
tral weight located for k ^ kp? The surprising answer 
is that the incoherent part of the u vs. k dispersion goes 
like — e(k) = — fc^/2m, "bending back" away from the 
chemical potential at large k. While this is expected in 
BCS theory and its generalizations for a paired super- 
fluid, we argue that this unusual dispersion is a universal 
feature of all dilute Fermi gases, even those with a normal 
(non-superfluid) ground state. We find that the spectral 



weight of C/fc in A(k, w) is centered about i 



-e(k) in 



a range of energies of order vpk for normal Fermi gases. 
Most of the spectral weight (1 — C/fc^) is, of course, cen- 
tered about w « +e(k), but these states are not occupied 
and do not contribute to n(k). 

This bending back is clearly visible in the data of 
ref. [H for attractive fermions near unitarity and near 



or above Tc- However, it is hard to separate the effects of 
the finite temperature pairing pseudogap and normal 
state interaction effects. In particular, a bending back of 
the dispersion above Tc cannot by itself be used as evi- 
dence for a pairing pseudogap in view of the normal state 
results described below. 

We will first focus on two systems where the ground 
state is a normal Fermi liquid: (a) the hard-sphere di- 
lute Fermi gas, and (b) the highly imbalanced attractive 
Fermi gas. We then turn to the superfluid ground state, 
where we will argue that, in the BCS limit, the unusual 
dispersion is dominated by interaction effects rather than 
the effect of pairing. We conclude with implications for 
RF spectroscopy experiments. 

Dilute repulsive Fermi gas: Wc begin with the 
3 dimensional hard sphere Fermi gas with dispersion 
e(k) — k'^ /2m, mass m, density n = kp/iir'^, and scatter- 



ing length a > with na'^ 1. (We set h = ks = 1-) Its 
thermodynamic and Fermi-liquid properties were stud- 
ied by Galitskii and Lee, Yang and Huang; see Sec. 5 of 
The high- A: tail was also calculated [12: '^(k) — 
(kpa)^ (2/37r)^ (kp/k)'^. Here we compute ^(k,a;). 

In the low density limit na'^ <C 1, the most important 
physical process is repeated scattering in the particle- 
particle channel. The corresponding sum of ladder di- 
agrams r is given by T~^{Q) = 1/g — L{Q), where 
Q = (Q, iQi) with iQi = i2iTrT and the bare interaction 
g is related to a via 1/g = 77i/(47ra) — ^j^. l/[2e(k)]. Fur- 
ther L{Q) ^TY.^ G^{k + Q)G°{-k) where k = (k, ifc„) 
with ikn = i{2n+l)TTT and G°{k) = l/[iA;„ - ^(k)] is the 
bare Green's function with the energy ^(k) = e(k) — u, 
measured with respect to the chemical potential /x [14| . 
Note that one can obtain an analytically closed form ex- 
pression 15[ for L(Q, H. + iO~^). For the hard sphere gas 



we can make a further simplification T g + g^L. 



The Matsubara self-energy E(A;) : 
yields i;(k, ik^ ^ u; + iO+) = ReE 



rE,r(fc+g)GO(g) 
ilmE. where 



ims](k, Lo)^J2 i™r(Q, [e (-e(q)) - (i) 



at T = 0, with Q = k + q and n = lu + ^(q). ReS 
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FIG. 1: (color online) Logarithmic intensity plot of 
A{]i,uj)eF/{kFa)'^ for the repulsive Fermi gas [kpa — 0.1; 
na^ = 3.4 x 10^^). The most intense (red) line at tj « 
is the quasiparticle. We focus on the unusual dispersion 
centered about uj — — e(k) (black dashed line) in the range 
uj = — e(k) — Sep ± 2vpk (white dashed lines); see text. 



is obtained numerically by a Kramers-Kronig transform 
on ImE. The spectral function is then obtained us- 
ing A(k,uj) = — Im [w — ^(k) — I](k, w)] /tt and plot- 
ted in Fig. [T] on a log scale. We see the most intense 
feature, corresponding to the Landau quasiparticle near 
kp^ tracks u w up to many-body renormaliza- 

tions [l7| . However our main interest is in the much 
less intense, incoherent spectral feature that follows an 
Lu = — e(k) dispersion and dominates ri(k) at large k. 

To understand this "bending back", we write A sa 
|ImS(k, lli)|/ [7r(a; — e(k))^] . We need to determine when 
ImE(fc,a;) is non-zero for k^ kp and w < 0. To under- 
stand our result qualitatively, consider the diagram in 
Fig. O The dominant contribution comes from small val- 
ues of both IQI and fJ. (For large values of these variables 
there is no spectral weight ImF for two-particle scatter- 
ing.) Thus q ~ — k and a; ~ — C(q) — — for k >• kp. 




FIG. 2: (color online) Top: Kinematics of the processes that 
contribute to imaginary self-energy in eq. ImE is non- 

zero when the shaded rectangle (allowed by kinematics and 
thermal factors) overlaps with the region Q, > r2o(Q)(in which 
ImF is non-zero). This leads to the condition Qmin < Qo- 
Bottom: Diagram contributing to ImE. 



0.06 
0.05 
0.04 
0.03 
0.02 
0.01 




I 1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 






-(k/kp)\/(kpaf 






-- (2/331)^ 


1 




- 50 n,^ / (kp if 


3 


4 


5 6 7 8 ^ 
k/kp 



FIG. 3: (color online) Momentum distribution tail for the 
dilute repulsive Fermi gas with no? = 3.4 x 10^^. 



This shows that A ^ for uj around negative e(k). 

To make this more quantitative, we use eq. ([1]). From 
the structure of L(Q,fi), it follows that Imr(Q,fi) ^ 
when n > QoiQ) = minp{^(p + Q/2) + ^(-p + Q/2)} = 
e(Q)/2 - 2fi; see Fig. H From the difference of 9- 

1/2 

functions, kp < q < (7max('^) = kp[l + \uj\/€p] . This 
implies that — |w| < < 0. Together with the kine- 
matical constraint \k — q\ < Q < k + q, this leads to 
Qmin = <?max(w)|. For non-zero ImE we thus need the 
kinematically allowed region (shaded rectangle in Fig. [2]) 
to overlap with > ilo{Q). This leads to the simple con- 
dition Qniin < Qo, where the definition fl{Qo) ~ leads 
to Qo = 2k p. (We have also found, but do not discuss 
here, the w > threshold for Ay^ 0.) 

Solving \k — qmax{^)\ ~ "^kp, we find ^(k, a; < 0) 7^ 
in the range of energies uj = ~e(k) — Sep ± 2vpk; see 
Fig. [H For k ^ kp, this simplifies to |a; + e(k)| < 2vpk. 
Although the width of this range grows linearly with fc, it 
becomes small relative to the central energy which grows 
like —k^ for large k. We plot in Fig. [3] the n(k)-tail 
using dujA(k, oj) and find that it agrees with the 
analytical result [l2|. The incoherent spectral weight in 
A{k ^ kp,uj) in the interval |a; -|- e(k)| < 2vpk is thus 
precisely {2kFa/ ?>t:)'^ {kp / k)'^ . 

Highly Imbalanced Fermi gas: We next turn to 
a two-component attractive Fermi gas with scattering 
length a tuned through a broad Feshbach resonance [l^ . 
While the ground state for equal spin populations is a su- 
perfluid exhibiting the BCS-BEC crossover, we consider 
the different regime of large spin imbalance n^/n^. There 
is by now considerable theoretical [T9l.[20| and experimen- 
tal evidence 0, [1| that, for sufficiently large imbalance, 
superfluidity is destroyed for a large range of values of a, 
including unitarity \a\ = 00, and the ground state is a 
(partially polarized) normal Landau Fermi liquid. 

For large |a|, we use the number of fcrmion species 
2M with an S'p(2A/')-invariant interaction as an artificial 
parameter to control the calculation in a large- A/" expan- 
sion [1^, [i^l- To first order in l/A/", ladder diagrams in 
the p-p channel determine the self-energy. The resulting 
expressions arc similar to those used above and we show 
them schematically, highlighting the differences that arise 
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FIG. 4: (color online) Logarithm intensity plot of Ak{uj)epf 
for minority particles in the unitary Fermi gas with imbalance 
ni/n-f — 0.01. The white dashed lines cj = — e(k) iawFt^ S'l'e 
derived in the text; the black dashed line is a; = — e(k). 



from spin-imbalance. We now have L ~ TJ2 G'ijJGj where 
Glik) = " CaM with ^,(k) - e(k) - The 

minority self-energy is given by Sj. = FG^. ImSj^ is 
then given by eq. ([T]) with ^ replaced by both in the 
0-function and in the definition of il. 

We can analytically determine the energy range for 
which ImE, and hence A, are non-zero. In Fig. [2] we 
must now use fiof^ = e(Q)/2 — 2^ with 2/i = fif + fi^. 
The final result [21| is that, for k ^ kp and w < 0, 
A(k, uj) can be non-zero only in the range of energies 
\io + e(k)| < avpfk where a 

CO. 



For concreteness, we focus here on unitarity \a\ 
y4(k, uj) for the highly imbalanced {ni/n-^ = 0.01) unitary 
gas is shown in Fig. [4l We have also verified that we 
get a 1/fc^ tail for n(k) in this system. Our calculation 
of yl(k, w) is controlled only within the l/A/'-expansion. 
Note, however, that the singularity structure in the large 
k limit is determined only by short-distance properties of 
the two-body problem in vacuum, while the strength of 
the singularity C depends on the many-body state. The 
ladder approximation is exact for the two-body problem. 
Thus we expect the bending back in the spectral function 
to be robust beyond the l/A/'-expansion. 

Superfluid State: We now turn to a discussion of the 
superfluid ground state for a system with equal densities 
of up and down spins and an interaction described by 
a scattering length a. Unlike the normal Fermi liquids 
described above, a branch of the dispersion that tracks 
— ek at large k is very natural for the fermionic excitations 
in a superfluid [l^l- Nevertheless, even in this case, our 
analysis gives important quantitative insights. 

In BCS mean field theory the spectral function 
Amf{K^) = vlS{uj + E{k)) + ulS{u} - Ei}^)) where 
«k = 1 - = [1 - C(k)/i^(k)]/2. The excitation en- 
ergy £'(k) = ■\/^2(k) -I- with A the energy gap. For 
k > kp, Ei^) w e(k) and « AV2e^(k), so that 
AMF(fc » /cf,w < 0) R:! [AV2e2(k)](5(tj + e(k)). Thus 
we see that particle-hole mixing in the superfluid ground 



state naturally leads to a bending back of the dispersion. 

However, there is a (large) quantitative problem with 
this result even in the BCS limit {Vjkpa <^ —1), where 
one might have expected it to be the most accurate. Us- 
ing n(k) = (iwA(k, w), or directly from BCS theory, 
one finds that the momentum distribution nMF(k) — 
vl « AV2e2(k) = Cup/k^ for k > kp. The prob- 
lem is that the contact estimated from BCS theory 
Cmf ^ A^ ~ cxp{'~l/ kp\a\) is exponentially small in 
\a\. However, the exact answer [1, [l3l in the BCS limit 



is G = 47r 



as a — > . To understand why BCS 



theory gets the wrong answer for G we use the adiabatic 
relation Q d£/da = C / {^-Kma^) . As shown in [23j . 
interaction effects lead to power-law corrections in |a| in 
the ground state energy density f , which are numerically 
much more important than the essentially singular cor- 
rections coming from pairing. In the extreme BCS limit, 
the contact is dominated by the Hartree term in £ with 
calculable corrections (23| . 

Thus the actual A{k kp,uj < 0), even in the BCS 
limit, is dominated by interaction effects beyond BCS 
mean field theory. This results in a spectral weight G ~ 
|ap arising from interaction effects which exist even in the 
normal state, rather than resulting from pairing, which 
only makes an exponentially small contribution. 

Implications for RF spectroscopy: The physical 
effects we have discussed above lead to directly observ- 
able consequences in RF spectroscopy experiments where 
an RF pulse is used to transfer atoms from one hyper- 
fine level to another. The interpretation of these experi- 
ments is often complicated by two difficulties: the inho- 
mogcneity of trapped gases and severe final state inter- 
actions. The first problem has been solved in the usual 
("angle-integrated") RF experiments using tomographic 
techniques. Final state effects are not an issue in 
0, and have been eliminated in ^Li by suitable choice 
of hyperfine levels 0, Q • We emphasize that but for this 
it would be very difficult to disentangle strong interac- 
tions in the many-body state (self energy effects) from 
final state effects (vertex corrections) [2J]. We will thus 
work in the (now experimentally relevant) limit where we 
ignore all final state interactions. 

Linear response theory then leads to the RF absorption 
intensity /cr(k,a;) = A<j(k,5<^(k) - a;)/(^cr(k) - to) where 
Lo is the RF shift. The Fermi function /(w) ensures that 
only occupied states can be excited by the probe. We 
set the multiplicative factor of the RF matrix element to 
unity so that J duj Ia{}^, ^) = 

Angle-resolved RF experiments [sj directly probe (the 
occupied part of) the spectral function A(k, ui) and can 
see its unusual dispersion for fc fcj?, w < 0. As already 
noted in the Introduction, one cannot identify this bend- 



ing back with pairing pseudogap physics |ll| since this 
universal feature also occurs in normal Fermi liquids. 

The consequences of our results for the angle-averaged 
RF experiments, which measure Ia{^) ~ -^f^C^' '^)' 
are more subtle. We now show that the unusual disper- 
sion at large k has a quantitative effect on the prefac- 
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tor of the universal high-w tail [25| in Ig-. We rewrite 

lA^) = Ek I-oo d^A^iK n)S{fl - e<.(k) + w) at T = 0. 
In the w — > oo limit, large negative fl values, centered 
about r2 « — dominate. We thus find /o-(w — > oo) ~ 
X;k«^(k)<5(a;-2ek). Using n^(k) « C/fc* for fc > fc^. 
we thus find that jH] /<,(w ^ cxd) « (C/47r2VT^)a;-3/2^ 
The characteristic power-law is independent of the phase 
(normal or superfluid) of the Fermi gas, though the value 
of C does depend on the many-body state. 

The a;"'^/^ tail has been discussed by various authors 
[12: m, [131 . First, we have derived this result analytically 
under very general conditions [1^. Second, our results 
show that one must be very careful in interpreting its 
origin. We emphasize that this tail and the bending back 
of the dispersion, to which it is intimately related, arise 
from short-distance "contact" physics, and should not by 
itself be taken as evidence for pairing. 

Conclusions: We have shown that there is an un- 
usual feature in the large momentum structure of the 
single-particle spectral function of all dilute Fermi gases, 
normal or superfluid, which is closely related to the uni- 



versal short distance features discussed by Tan and oth- 
ers [^ [l3] . This is an incoherent branch of the dispersion 
where w goes like negative ek [13 that is quite unex- 
pected in a normal Fermi gas. Nevertheless, this is ex- 
actly what we find in the two systems where the ground 
state is known to be a normal Landau Fermi liquid: the 
hard sphere Fermi gas and the highly imbalanced, attrac- 
tive Fermi gas. Even in a BCS superfluid, we show that 
this bending back at large k is dominated by interaction 
effects rather than the pairing gap. 

Note Added: Recently, we learned of a work by 
Combcscot. Alzetto and Leyronas [11], where the ap- 
proximation ImE(fc kF,u} < 0) cx i5(w -|- e(k)) is used 
which leads to a sharp feature in A(k, oj). While this 
may be sufficient for computing "integrated" quantities 
like n(k), it does not capture the incoherent structure in 
A(k, uj) described here. 

We acknowledge discussions with D. Jin and E. Taylor 
and support from Insf-dnir 07 06203 and ARO W911NF- 
08-1-0338. 
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